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Abstract 

^ ' We study the transversal spin-spin correlations in the non-equilibrium steady state of 

the XY chain constructed by coupling a finite cutout of the chain to the two infinite parts 
to its left and right acting as thermal reservoirs at different temperatures. We prove that the 
spatial decay of these correlations is at least exponentially fast. 
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1 Introduction 

The XY chain is the one-dimensional integrable spin system, introduced in [ 1 1 1 (see also J9)), 
whose formal Hamiltonian is specified by 

H =-\L{( 1 +r) <X! ( V +1) + (1 " 7) 4 X) 4 +1) +2A a 3 W } , (1.1) 

xeZ 

fx) 

where <7j denotes the Pauli matrix at site xeZin the transversal directions, j = 1,2, and in 
the longitudinal direction, j = 3. The parameter 7 G (—1,1) describes the anisotropy of the 
spin-spin coupling whereas A el stands for an external magnetic field. 

Since the discovery of their ideal thermal conductivity, quasi-one-dimensional Heisenberg spin- 
1/2 systems, made from different materials, have been intensively studied experimentally (see 
[fT4l IT31 : SrCu02 and Sr2Cu03 are often considered to be the best physical realizations of one- 
dimensional XYZ Heisenberg models). 

Not only this highly unusual transport property motivates the theoretical study of such non- 
equilibrium models, but the XY chain also represents one of the simplest non-trivial testing 
grounds for the development of general ideas in rigorous non-equilibrium theory. 

Already in ifTTl . for vanishing external magnetic field A = 0, the transversal spin-spin corre- 
lations in the ground state and at nonzero temperature could be expressed by means of deter- 
minants of large truncated Toeplitz matrices. Moreover, bounds on these correlations implied 
vanishing long-range order, at least in the isotropic case 7 = 0. Later, in lfT2l . this study was con- 
tinued yielding an asymptotic evaluation of the transversal correlations with the help of Szego's 
theorem. Afterwards, in [6], almost the complete phase diagram in 7 and A for the behavior 
of the correlation functions both in the longitudinal and in the transversal directions were ob- 
tained, and so for zero and nonzero temperature: The result of this study for the case of nonzero 
temperature is that all the correlation functions are asymptotically exponentially decaying with 
a decay rate which depends on the magnetic field A and the anisotropy 7. 

In this note, we study the large n behavior of the transversal spin-spin correlations, 

Cj(n) = co(ofof), 7 = 1,2, neN, (1.2) 

in the non-equilibrium steady state (NESS) CO constructed in [5] (and already in [3 1 for 7= 0) 
in a setting which has become to serve as paradigm in non-equilibrium statistical mechanics: 
a "small" system which is coupled to two infinite reservoirs which are in thermal equilibrium 
at different inverse temperatures /3l and fa. We prove that the spatial decay of Cj(n) is at 
least exponentially fast. This behavior contrasts with the one generally expected by the folklore 
which predicts quasi-long-range order out of equilibrium. 

In the following section, we give a brief informal description of our non-equilibrium setting 
for the XY chain. We refer to El|5l[l0| for a precise formulation within the framework of C* 
algebraic quantum statistical mechanics. 
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2 The non-equilibrium setting for the XY chain 

Consider the XY chain described by the Hamiltonian dl.ll) and remove the two bonds at the 
sites ±M, M > 0. Doing so, the initial chain divides into a compound of three noninteracting 
subsystems. This configuration is what we call the free system whose Hamiltonian 

Hq = H L + Hs + H R 

is built as in dl.ll) according to the union Z^ U Z$ U Zr, where Zi = {x E Z \x < — M — 2}, 
Z R = {xeZ\x>M+l},<mdZ s = {xeZ\-M<x<M-l}. The infinite pieces Z L and Z R 
will play the role of thermal reservoirs to which the finite system on Z$ is coupled by means of 
V = H — Hq. In contrast, the configuration described by H, i.e. the original XY chain on the 
whole of Z, is considered to be the perturbed system. 

In order to construct a NESS in the sense of lfT3*l . we choose the initial state C0q to be composed 
of (Xj, j8y)-KMS states C0j on Zj, j = L,R, and of the normalized trace state cos on Z5 (T/ denotes 
the time evolution generated by Hj, and /3, is the inverse temperature), i.e. we set 

It is well-known that the Jordan-Wigner transformation maps the XY model on a model of 
free fermions with annihilation and creation operators c x , c*. Using scattering theory on the 1- 
particle Hilbert space of these fermions, the NESS CO w.r.t. the initial state 0Jq and the perturbed 
time evolution t' have been constructed in [5|, 

oj(c x ) = lim (0D(T f (c x )). 

t— > +00 

It has been shown in [5] that CO is a quasi-free state with 2-point operator S, 

co(B*(f)B(g)) = (f,Sg), (2.3) 

where /, g G/ 2 (Z)®C 2 ~Z 2 (Z)©/ 2 (Z),and5(/)=I xeZ (/ + (x)c*+/_Wc,)for/= (/+,/_). 
Moreover, S has been explicitly computed in [ 5 1 . In the Fourier picture, / 2 (Z) ~ L 2 (T) (with 
T = {zeC||z| = l}), it reads 

S^)=(l+eP h ^ +5k ^y\ (2.4) 
with the parameters /3 and 8 given by 

P = l&R + h), S = (2.5) 
The 1 -particle operators h and k have the form 

h{%) = (cos^-A)®(7 3 -7sin^(g)a2, k(%) = sign(ir(£))/i(§) ® ct , 
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where the functions fc(^) and ju(^) are defined by 

Jc(^)=2Xsin^-(l-r 2 )sin2^ = ((cos^ - A) 2 + y 2 sin 2 ^) 1/2 . (2.6) 

Furthermore, Ob is the identity on C 2 , and o~i, 02, 03 denote the Pauli matrices, 





10" 




"01" 




" -/ ' 




" 1 




1 


, 01 = 


1 


, C7 2 = 


/ 


, ^3 = 


-1 



(2.7) 



3 Exponential decay of the correlation function 

Let us now state our result. 

Theorem Let CO be the quasi-free NESS characterized by (12.3b . and Cj(n), j = 1,2, the corre- 
lation functions (I1.2I) . For all values of the parameters < 8 < j8 < 00, 7 E (— 1 , 1), and A G R, 
?/ze transversal spin-spin correlations Cj(n) decay at least exponentially fast, and 

limsup g| J| < - / £ log [th(jSi/i(0/2)th(jSu/i(^)/2)] < 0. 
n— >°o n Z Jo In 

Remarks 1) If the temperature difference of the reservoirs vanishes, i.e. 5 = 0, the total system 
is in thermal equilibrium at inverse temperature fit = Pr, and our bound estimates from above 
the asymptotically exponentially decaying behavior from [6]. 

2) For 8^0, the correlation in the ^-direction decays like \/n 2 at infinity for all 7 G (—1,1), 
X G R, see f^. In contrast to this result, our theorem does not confirm the folklore about the 
change in the type of decay - from short range to long range - when passing from equilibrium 
to non-equilibrium. 

Before we begin the proof of the theorem, we rewrite the correlation function Cj(n) as the de- 
terminant of a In x In block Toeplitz matrix, in the same way as it has been done at zero and 
non-zero temperature in O HTJ . The reader not familiar with Toeplitz theory may consult Ap- 
pendix |X| before entering the proof. We restrict our attention to the transversal correlation func- 
tion C\{n) = C(n) in the 1 -direction (the 2-direction being analogous, see the Jordan-Wigner 
transformation for 02 in (13.81) below), 

C(n) = co(a[ 0) a { 1 " ) ), hgM. 

(x) (x) (x) 

The well-known Jordan-Wigner transformation expresses the spins o~j , , Oq , , x G Z, by 
means of fermionic creation and annihilation operators c| (= c*, c x ), 

a (x) = TS (x) + c *) j G (X) = . TS ( X ) ^ x _ ^ (x) = 2c * Cx _ h (3 g) 
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where = O3 . . . ojf ^ for x > 1, = 1, and = O3 . . . O3 for x < 1, see for example 
0. The element 7 1 stems from Araki's C* crossed product extension of the CAR algebra for the 
two-sided chain, see Q. It has the properties T 2 = 1, 7cq = — Cq7\ and Tel = cfj for x > 0. 
Plugging the Jordan-Wigner transformation (13.81) into the product ffj cj , we find 

where = c* ±c x . Moreover, we define r\± = (1,±1) G C and o4 ; = <5 X ® 7]± (with 5* G 
/ 2 (Z), 5jc(y) = 1 for y = x and zero otherwise). With the linear map B introduced in [1J, / 2 (Z) © 

Z 2 (Z) 3f= (/+,/_) h- £(/) = L xeZ (/+ (x)c* +/-(x)c,), we have = B{(^). Hence, we 
can express the correlation function C(n) as 

C(n) = (-l)' I a)(5(al V(4 1) ) 5 ( a - 1 V----^(«i^ 1) ) 5 (« ( ' ,_1) )5(a? ) ))- 
Since (0 is quasi-free |5|, the correlation function C(«) is a Pfaffian, 

C(n)=PfQ(n), 

where the matrix (Q.(n) jk)j n k=l G C 2nx2n is defined to be skew-symmetric, and, for j < k, 

£l{n) jk = (o{B{f } )B{f k )) with f 2i _ l = at 1 \f2i = cc i i\ i = l,2,..,n. 

The Pfaffian PfA of a matrix A G C 2 ^ 2 " is defined by PfA = I^sign(^)nLl A ^-i,^ with 
the sum running over all n in the permutation group £2,, which satisfy %2 kl fl^t+i > ?hk-l- If -A 
is skew-symmetric, A T = —A (A r denotes the transpose of A), the Pfaffian of A is related to the 
determinant of A through (PfA) 2 = detA. Thus, we are led to study the large n asymptotics of 
the determinant of Q.(n), 

C(n) 2 = deta(n). (3.9) 

In order to cast Q.(n) into the form of a truncated block Toeplitz matrix, we compute the matrix 
entries 

a£ ± = a>(B(a^)B(a^)), \ <j<k<n. (3.10) 
For this purpose, we rewrite the 2-point operator S(£ ) on L 2 (T) ® C 2 from (12.41) as 

3 

using the Pauli matrices from (12.71) . The first component sq(£) looks like 

= i + irignC^))^), (3-12) 
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and ,5(1) = (ji(£),J2(£),J3(£)) has the form 

s(£) = iVpiZMZ), KS) = ^(0,-ysin^cos^-l), (3.13) 

with the functions jc(^) and ju(^) from (12.6b . Moreover, we used the definition 

sh(a/i(£)) 



ch( J 8 i u(^))+ch(5 i u(^)) : 



Let us start with the computation of in (I3.10I) . Using the property B*(f) = B(Jf) for 
J : f = (/ + ,/_) i — > (/-,/+) from [5|, and the expressions in (13.111) . (13.121) . and (13.131) . we find, 
for 1 < j < k < n, and with the definition ej(^) = e~'^, that 

A++ = (ej®T} +1 Se k ®T} + ) L 2^ C 2 
27T 

The computation of A^ k _ yields A^ k _ = —A+ + . Next, the entry A+_ looks like 

A +- = {ej®T)+,Sek®?)-)i2®c 2 



a e 



(3.14) 



2 /o gte(£)+M£)K 

2jr <i^ cos<^ — A — /ysin^ 



wit: 



<k-M 



Finally, A J ^ + = —A+_ . Due to translation invariance, we can define A± ± = A J 2± X with x G Z. 
Using Ai_ = — A^ + , Ai + = — A~£_, and the symmetry property A+i = — A^ + , we can write 
in (13.91 ) in the form of a 2n x In truncated block Toeplitz matrix with 2x2 blocks a x , 



Q(n) 



where the blocks are given by 



a i ciq 

U-n— 1 ^n— 2 



Ai 



-A 



-Ai 



fl -(«-i) 

«-(n-2) 
«0 



, x G Z. 



(3.15) 



(3.16) 
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Using the relations a x = j Q l7t dt, / (2n) a(£,)e from (IA.20I) . we can extract from (13.161) the 
symbol a(£) 



sign(jc(0)95(0 -9(09)8(0 
9(0^(0 -sign( K (0)%(0 



(3.17) 



Here, <p a (<^) is defined in (13.141) . ?c(0> M(0 m A2.6I) . and #(£) is the complex conjugate of 

cosg-A-iysing - g 
^ M(0 

Therefore, denoting by T[a] the block Toeplitz matrix associated with the symbol a given by 
(13.171) . we can write Q.(n) as a truncated block Toeplitz matrix (see dA.231 )). 



Q(n) = T n [a]=P n r[a]P n . 



Proof of the Theorem For n e N, let sy < < ■■■< s^} be the singular values of T n [a], and 
denote by N en the number of singular values in [0, e], < £ < 1, including multiplicity. From 
the above computations, we have 



In In 

\C(n) | 2 = | detCl(n) | = | det T n [a] | = f[ sj < £^>" [ ] s) 



Thus, 

2n . . 2n 



io g |detr n [ fl ]| < £ io g 4 M) <£^(4 n) )iog4 n) 

j=N e .+i j=i 



where < ^ e < 1 is a smooth characteristic function with support in (e, \\T[a] \\ + 1) and value 
1 on [e + £ 2 , ||r[a]||]. Since the symbol a of T„[a] has the two singular values <pp — (p§ and 
q)p + cp§, the Avram-Parter theorem JH1 p.203] (see also (IA.26I) ') yields 

limsup-log|detr B [a]| < / -^tr[(# e log)(|a(OI)] 

n-^oo ft JO 27T 

= ^ 2 ^U £ log)(^(0- 95(0) + Ofclog)(9p(0 + 95(0)- 

Since the inequality holds for all £ > 0, and since for all £ E [0,2n] we have < <p«(0 — 
95(0 < 9/3 (0 + 95(0 < 1. we obtain 

limsup-log|detr„[a]| < / log [th(j3 L /i(0/2)th(j3*/i(0/2)] (3.18) 

n-s-oo n JO m 
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which concludes the proof. 



□ 



Remarks 3) We do not provide a lower bound in our proof (but see remark 1 for 5=0). Due 
to the lack of results on the invertibility of block Toeplitz matrices T[a], the behavior of the 
singular values ofT n [a] close to zero is, in general, rather difficult to control. 

4) In the critical regime of the XY chain, 7=0, A G [—1,1], and y ^ 0, A = ±1, where the 
nonnegative function jii(^) from (12.61) can have zeroes on [0,2tt], the bound (13.181) remains 
finite due to the integrability of \ogx at the origin. 

5) A straightforward proof using the inequalities | det T n [a] \ < \\T n [a] \\ 2n and 

H^nHH < II^HII = ess sup ||tf(£)||_^(£W) = ess sup S2[a{£>)) 

(see and (IA.24I) - (IA.25I) J yields the weaker bound \ detT n [a] \ < th 2 "(j3/?||jii||oo/2) which, in 
contrast to (13 .18b . is valid for all n. On the other hand, in the limit [5r — > °°, this bound is void 
whereas (13.181) still implies exponential decay as long as < °°. 



A Toeplitz operators 

Toeplitz operators ||8l p. 185] Let N EN. We define the space 1% of all C^-valued sequences 
/ ={fi)U, C N , by 



1/2 



# = {/:N-C 



N i 



<oo}, 



(A.19) 



where || • || c at denotes the I norm on C*. We write I 2 = If. Let {a x } xeZ be a sequence ofNxN 



■*NxN 



matrices, a x G The Toeplitz operator defined through its action on elements of l^ by 

/ l— *■ {E7=i a i-jfj}7=\ 1S a bounded operator on /jy, if and only if 

^ d$ 



for some a G L^ xN (see p. 186]), where we define (with T = {z G C | \z\ = 1}) 



^NxN 



{0:T 



iNxN 



(j) i jeL 00 (T),i,j=l,...,N}. 



(A.20) 



(A.21) 



In this case, we write the Toeplitz operator as 

T[a) = 



uq a— i a— 2 
a\ ag a— i 
a2 a\ ag 



(A.22) 
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The function a G L^ xN is called the symbol of T[d\. \f N = I the symbol a G L°° = L^ xl and 
the Toeplitz operator T [a] are called scalar, whereas for N > I they are called block. 

For nGN, let P n be the projections on ijj, 

With the help of these P n , we define the truncated Nn x Nn Toeplitz matrices as 

T n [a]=P n T[a]P n \ lmPn . (A.23) 

Norm of Toeplitz operators [ 8 , p. 1 86] The norm of a Toeplitz operator is related to its symbol, 

||r[fl]|| = ||a||o., (A.24) 

where ||a||oo is defined to be the operator norm of the multiplication operator acting on L 2 (T) 
by multiplication with the matrix function a G L^ xN . From p. 95], we have 

||a||oo = esssup ||a(£)|| (A. 25) 

where || • || _g(vN\ is the operator norm induced by the Z 2 norm on C . 

Avram-Parter theorem [8, p. 203] Let s\ , be the singular values of the truncated 
Toeplitz operator T n [a] with block symbol a G L% x m, and g a continuous function with com- 
pact support, g G Cq(IR). Then, 



1 Nn i \ 1 r l7Z dt 
lim — Y t g{sf ) ) = ± ^trg(K^)|), (A.26) 



where trA is the trace of the matrix A G C mxm , m G N, and \A\ = yA*A its modulus. 
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